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ABSTRACT 

We extend the mass-halo formalism for analytically generating power spectra to al- 
low for the different clustering behavior observed in galaxy sub-populations. Although 
applicable to other separations, we concentrate our methods on a simple separation by 
rest-frame color into "red" and "blue" sub-populations through modifications to the 
(N)(M) relations and halo distribution functions for each of the sub-populations. This 
sort of separation is within the capabilities of the current generations of simulations as 
well as galaxy surveys, suggesting a potentially powerful observational constraint for 
current and future simulations. In anticipation of this, we demonstrate the sensitivity 
of the resulting power spectra to the choice of model parameters. 
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1 INTRODUCTION 

The difference in clustering between intrinsically red and 
blue galaxies has been known since galaxies were first clas- 
sified into types by Hubble. More recently, the state of cos- 
mological simulations has reached the level where galaxy 
evolution can be modeled semi-analytically (White & Rees, 
1978; White & Frenk, 1991) to produce simulated catalogs 
with realistic distributions of galaxy colors and types (Kauff- 
mann et al., 1999; Somerville & Primack, 1999; Benson et al., 
2000;). Within the framework of these types of simulations, 
the power spectra and bias of "red" and "blue" galaxies have 
been measured and found to be in relatively good agreement 
with previous data measurements. In this paper we seek to 
derive an analytic method for generating power spectra and 
biases for these red and blue galaxies from within the mass- 
halo model. 

The first papers developing the revisions of the orig- 
inal mass-halo formalism for dark matter and galaxies fo- 
cused on real-space correlation functions (Sheth & Jain, 
1997; Jing et al., 1998; Peacock & Smith 2000). Recently, 
these treatments have been extended for the calculation of 
power spectra (Seljak, 2000; Scoccimarro et al., 2000; Ma & 
Fry, 2000). Although the results generated by these models 
closely matched the observed power law behavior for galaxy 
power spectra (Hamilton et al., 2001), they did not address 
the observed different clustering within sub-populations. We 
build on the results from these earlier power spectra treat- 
ments to produce physically motivated models for power 
spectra and biases for galaxy sub-populations. In an ef- 
fort to keep our model as general as possible, we will re- 



frain from specifying precisely what constitutes a "red" or 
"blue" galaxy; our approach should work equally well for 
galaxies separated by type or any other observable with dif- 
ferent clustering properties. 

In we review the formalism for calculating galaxy 
power spectra. With this laid out, §[] details the modifica- 
tions to the halo profiles and the number-mass relation for 
galaxies populating the halos. Additionally, we work through 
the modifications to the formulae from ^ necessary to calcu- 
late the sub-population power spectra and cross-power spec- 
trum, as well as the associated relative biases. Finally, in 
we explore the power spectrum space spanned by the new 
parameters added to the standard model in §H. 



2 CALCULATING POWER SPECTRA 

Following the treatment in Seljak, we have four essential 
ingredients in our galaxy power spectrum: the halo profile, 
the halo mass function, the halo biasing function, and the 
galaxy number function. Once these four components are 
determined, either from observations or simulations, we can 
fold them together to produce the predicted power spectrum 
for that particular model. 

We begin by generating a halo profile, p(r), which 
is parametrized along the lines of the profile derived by 
Navarro, Frenk and White (1996; NFW, hereafter), 



p(r) = 



PS 



(r/r s y a (l + r/r s ) (3+t *> 



(1) 



where rs is the universal scale radius and ps = p( r s)- In 
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most halo-model treatments, rs is replaced by a concen- 
tration, c = r v /rs, where r v is the virial radius defining 
the region where the fractional overdensity of the halo S v 
is approximately 200 and c is a weak function of halo mass 

(c = c (M/M*f, where c ~ O(10) and (5 0(10^)). 

The traditional NFW profile gives a = — 1, while the Moore 
profile (Moore et al. 1998) has a = -3/2. We will use 
a = —1.3 for the calculations in this paper, but the gen- 
eral results are largely insensitive to the choice of a. Bullock 
(2001) gives Co = 9 for a pure NFW profile; using Peacock 
& Smith's relation, Co w 4.5 for a Moore profile. Since we 
are using an intermediate value of a, we choose Co = 6 and 
(3 — —0.15 for all the calculations in this paper. In principle, 
one can also consider scatter in the concentration at a given 
mass (Scoccimarro, et al., 2001), leading to an integral over 
the distribution of c, but the magnitude of this effect is small 
enough that it can be safely ignored. 

Rewriting Equation [l] in terms of the concentration and 
the mass, we get 



p(r,M) = 



(rc/r v ) Q (l+rc/r„) 



3+a ■ 



where 



Ps 



3M 
4ir8 v p' 

S v pc 3 (M) 



:(M) 



d X 



(i + x) 



3 + a 



(2) 



(3) 



(4) 



p is the mean matter density and M is the mass of the halo. 
Since we will be working in wavenumber space when we 
generate the power spectrum, we actually need to consider 
the Fourier transform of the halo profile, 



y(k,M) = 



4nr ■ p(r,M) 



sin(fcr) 
kr 



dr. 



(5) 



where we have normalized over mass so that y(0, M) — 1 and 
y(k > 0, M) < 1. Note that this implies that p(r > r v ) = 0, 
truncating the mass integration at the virial radius. 

With this in hand, we can move on to the next compo- 
nent of the halo model, the halo mass function, (dn/dM). 
Traditionally, this mass function is expressed in terms of a 
function f(v), 

4^dM = JLUv)dv, (6) 
dM M jy ' ' w 

where v relates the minimum spherical over-density that has 
collapsed at a given redshift (S c (z), 8 C (0) — 1.676 for an 
f^m = 0.3, S1a = 0.7 cosmology) and the rms spherical fluc- 
tuations containing mass M (a(M,z)) as 



S c (z) 
<r(M) 



(7) 



We define M* as the mass corresponding to v = 1. The 
functional form for f(u) is traditionally given by the Press- 
Schechter function (1974). This form tends to over-predict 
the number of halos below M„, so we use the form found 
from simulations by Sheth and Torman (1997), 



v.f{v)~{l + v'- p ) 



t-ps ,1/2 - 

v )u e 



//2 



(8) 



where v 1 = av, a — 0.707 and p = 0.3. This relation is 
normalized by requiring that 



dfi 

j^MdM = I f{v)dv = 1 



(9) 



for the dark matter distribution. 

On nonlinear scales, we expect the halos to cluster more 
strongly than the mass, and vice versa for linear scales (Mo 
& White, 1996). This means we need to positively bias the 
clustering of the high mass halos relative to the low mass 
halos. We can generate this sort of halo biasing scheme for 
the ST mass function using 



b(u) = 1 + V — 1 + 



2p 



(10) 



In order for the eventual power spectrum to reduce to a 
linear power spectrum on large scales, we need to impose 
the further constraint that 



f(v)b{u)dv = l, 



(11) 



requiring that the biased halos with mass greater than M* 
be balanced out by anti-biased halos with mass less than M, . 
This integral is satisfied automatically if we use Equation |Io| 
and have properly normalized f(v). 

Using just these three components, we can generate the 
power spectrum for the dark matter. However, in order to 
predict the galaxy power spectrum, we need to know how 
many galaxies are in a given halo (under the assumption that 
the distribution of galaxies in the halo follows the halo pro- 
file). Currently no theory completely informs the formation 
of galaxies given a halo mass, but the traditional form of the 
(N)(M) relation (Jing et al., 1998; Kauffmann et al., 1999; 
Benson et al., 2000; White et al., 2001) has the galaxies pop- 
ulating the halo as a simple power law, (N)(M) ~ (M/Mq) 1 , 
where Mo sets the unit mass scale and 7 < 1. Additionally, 
one can put in constraints on the minimum mass to form a 
galaxy and other modifications. For the purposes of the for- 
malism for calculating the power spectrum, however, we can 
put aside the question of precisely what this function looks 
like. The inclusion of galaxies does change the normalization 
of Equation |^ to 



——f(v)dv = -, 
M K V ) P 



(12) 



where n is the mean number of galaxies. 

On large scales, the power spectrum is dominated by 
correlations between galaxies in separate halos. We need to 
convolve the halo profile with the mass function to account 
for the fact that halos are not pointlike objects. Since we 
are in Fourier space, we can perform the convolution using 
simple multiplication. The halo-halo power {Pgg{^)) i s then 
simply, 



Ptcik) = PLIN(fe) 



(N) 
M(u) 



b(v)y(k,M)dv 



,(13) 



where PlinW is the linear dark matter power spectrum. In 
the small k limit, this reduces to a simple linear bias ((b)), 



(14) 



For small scales, the dominant contribution to the 
power spectrum comes from correlations between galaxies 
within the same halo. This single halo term is independent 
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of k at larger scales, giving it a Poisson-like behavior. In or- 
der to account for the fact that a single galaxy within a halo 
does not correlate with itself, we use the second moment of 
the galaxy number relation, ((N(N — 1))) to calculate the 
Poisson power (Pgg(^))j 

p ° c(fc) = r m m ^ i)) uk ' m<dv ' (i5) 

Seljak takes ( = 2 for (N(N - 1)} > 1 and ( = 1 for 
(N(N — 1)} < 1; this is done to account for the galaxy at 
the center of the halo in the limit of small number of galax- 
ies. In the limit of large numbers of galaxies, (N(N — 1)) 
approaches (N) 2 , but in the small number limit, (N(N — 1)) 
can be approximated by a binomial distribution. This can 
be implemented by using the fits given in Scoccimarro et al. 
(2000), letting 

(N(N-l)) =a 2 M (N> 2 , (16) 

where a M = 1 for M > lO 13 h~ 1 Af and a M = 
log( A /MVlO 11 M ) for M < 1O 13 /i _1 M . Clearly, this 
will not be exact for an arbitrary (N)(M), but it should be 
close enough for our purposes. Adding Po%{k) anc ^ ^ggWj 
we recover the galaxy power spectrum at all wavenumbers, 
P GG (k). 



tively) would be the simple power laws alluded to earlier, 

In practice, this form actually over-determines the functional 
form of the power laws; since we know the functions pass 
through unity at some mass scale, we can set JVro = -/Vbo = 
1 and determine the relative contribution of red and blue 
galaxies by selecting Mro and A/bo appropriately. 

Such a model does reasonably well, but the GIF simula- 
tions (Kauffmann, et al., 1999) point to an extra abundance 
of blue galaxies at small halo masses (M ~ 1O 12 /i _1 M ). 
Sheth et al. (2001) include this effect by the addition of a 
Gaussian term to the (N)b relation, 

(N)b(M) = (^) 7B +Ae-^°^- M ^ 2 , (18) 

where A is 0(1/2) and Mb s is the logarithm of the mass 
corresponding to the peak in the Gaussian component. All 
told, this gives us six tunable parameters for the (N) rela- 
tions. For the purposes of our power spectrum calculations, 
we will add an additional lower mass cut-off for the (N)r re- 
lations at lO 11 /t~ 1 M ; the resulting power spectra are not 
terribly sensitive to the precise value of this cut-off. 



3 GENERATING RED &: BLUE POWER 
SPECTRA 

Within the framework presented above, there are a number 
of parameters which might be modified to generate mod- 
els of different power spectra for red and blue galaxies. One 
could modify the concentration index for each galaxy pop- 
ulation, change the halo biasing relation, etc. In this paper, 
we focus on two modifications to the standard treatment: 
a modification of the (N) relations and the halo profiles for 
each of the galaxy types. 

The physical motivation in both of these cases is clear. 
Semi-analytic models for galaxy formation indicate that the 
primary determinant of galaxy color is the epoch of ini- 
tial gas cooling to form the initial stellar population. Red 
galaxies tend to form earlier, appearing in the deepest over- 
densities, while the current blue galaxies form later when gas 
in the shallower potentials and outskirts of the larger poten- 
tials has cooled. Given this difference in development, the 
prospect that the efficiency of galactic formation (and hence 
number of galaxies produced within a halo of a given local 
halo mass) would vary for the two epochs is a reasonable 
conclusion. Likewise, for the different halo profiles, we know 
from observations as well as simulations that red galaxies 
tend to populate the centers of galaxy clusters and filaments, 
while blue galaxies are more numerous at the fringes of struc- 
ture and in the field. Changing the distribution function for 
the different colored galaxies within a halo to reflect these 
observations is an obvious step. 

3.1 Modifying (N)(M) 

The simplest model we can adopt for the galaxy number re- 
lations for the red and blue galaxies ((N)r and (N)b, respec- 



3.2 Modifying Halo Profiles 

In modifying the distribution functions for red and blue 
galaxies, we have a number of constraints. First, we require 
that the sum of the matter associated with red galaxies and 
that with blue galaxies match the total distribution of mat- 
ter at all halo radii. This is not to say that all the matter 
ends up in galaxies; rather, it merely requires that the com- 
bined distribution of red and blue galaxies in the halo match 
the total galaxy distribution. Although this would appear to 
suppress natural correlations between red and blue galaxies, 
Sheth & Lemson (1999) show that this sort of clustering by 
conservation of number works reasonably well. 

There are any number of profiles we might choose to 
consider for the red and blue sub-profiles. In principle, the 
shapes of these profiles could be found from analyzing the 
results of simulations, but for the purposes of this exercise 
we will forgo that complication. Instead, we will restrict our- 
selves to profiles of a similar form as that given in Equation [j] 
with different values of ps and a for the red and blue sub- 
populations. Of course, the sum of two profiles with differing 
values of a do not quite match a profile with a third value 
of a, but we can come reasonably close with clever choices 
of «r and qb for the red and blue populations, respectively. 

In order to find the proper values for our extended halo 
parameter set, we need only know the relative abundance of 
red and blue galaxies at two points along the radial profile. 
At large radii relative to the virial radius, all of the profiles 
go as p ~ r~3, meaning that the ratio of the number of blue 
galaxies to red galaxies (rj) in this regime gives us the ra- 
tio of psb/psr immediately, once we have transformed from 
number of galaxies to mass assigned to galaxies (see below) . 
We consider the halos truncated at these large radii when 
we calculate the power spectra, but this allows us to set the 
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relative normalization of psr and psb under the assumption 
that 7] does not vary much beyond r v . 

With this ratio determined, the choice of the second 
radial point is somewhat arbitrary, but clearly we would 
prefer this point to be at very small radius where the profiles 
go as p ~ r a . For all of the calculations, we set the inner 
radius such that ric/r v = 0.1. This demands a resolution 
smaller than most simulations can provide at the moment, 
but does guarantee that we are well into the p ~ r Q regime. 

For a given ratio (p) of red to blue galaxies at small 
radius, the difference in a for the red and blue profiles is 



Act 



«R — 



log(l + nc/r v ) - log(nc/r v ) : 



(19) 



where we have assumed that the mass follows the number 
of galaxies linearly. However, to properly conserve the mass 
in the halo, we need to transform p and r\ into mass-space 
using the (N) functions for the red and blue galaxies to yield 
p! and rj' , respectively. In the absence of the Gaussian term 
in the (N)b relation, these transform as 



V 



'/ 



i/ 7B Mbo 
M R o 
i/ 7R Mrq 
Mbo' 



(20) 



If we wish to include the Gaussian term, a bit of algebra 
gives the slightly more complex form of the transformation: 



V = 

I 

P = 



r) — Ae 



4(log(M R0 )-Af B3 ) 



i/7b M B 



M R0 

Mrq 

Mbo' 



(21) 



As one might expect, the limit that Mro,Mbo 2> 10 Mbh , 
Equation |il] effectively reduces to Equation Once we 
have made this transformation, we then replace p and r\ 
in Equation [l9] with p! and rj . With this relation between 
qr and ob in hand, we can perform a simple search over 
values of qr to find the sub-profiles that combine to closest 
match an overall profile with a given value of a. Since we 
know that Aa must be positive, this relation guarantees a 
flatter distribution of blue galaxies in the center of halos 
relative to red galaxies, something which would be much 
more difficult to accomplish had we varied c for the differ 
galaxy populations. 

Using psr, Psb, Qr and a B , where 



PSR 



Psb = 



-ps 



(22) 



V 



TPS, 



rj + 1 ' 

we can recalculate y(k,M) along the lines of Equation ^| 
making i/r and j/b- We can also recover the distribution of 
galaxies by color/type within a given halo similar to that 
seen in simulations (cf. Diaferio et al., 1999), as seen in Fig- 
ure Clearly, the agreement is not perfect, but it is close 
enough that we can proceed with confidence that the basic 
approach is reasonable. These diagrams inform our fiducial 
choices of p and rj such that p ~ O(10) and rj ~ O(a/(10)); 
for our models we take p = 10 and rj — 4. 
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Figure 1. Fraction of "red" and "blue" galaxies as a function of 
rc/r v compared to the distribution of galaxies from the GIF sim- 
ulations. For the simulated data, red galaxies were taken to have 
g' — r' < 0.6 in the rest frame for SDSS magnitudes (Fukugita, 
et al. 1996). 



3.3 Calculating Subpopulation Power Spectra and 
Biases 

Plugging all this into the standard formalism, we can calcu- 
late the power spectra for these galaxy populations, as well 
as their cross-correlation, in the same manner as described 
in Seljak. In order to calculate the power spectra properly, 
we need to re-normalize f(v) for each sub-population using 
Equation [l2| to account for the differences in (N) and n: 



(N)r 
M{v) 

(N)b 
M(u) 



fa.(v)dv 



rtR 
P 

P 



(23) 



Once this has been done, we can insert the above (along with 
the color- dependent halo profiles) into Equations |l3| and |l5| 
to generate the power spectra for red and blue galaxies: 



_ f 

PLm(k) 


_£_ 

fin 


PrrW = P 


( 


PBB(k) \ 


P_ 

flB 


Pisik) =p 


p 00 

Jo 



<N)r 



b{u)y n {k,M)dv 



\y K {k,M)\ c dv, 



y ' M{v 
(N(N-1)) R 

lB{v)^^b(y) V B(k,M)dv 
(N(N-1)> B 



(24) 



(2n) 3 n 2 B M(u) 



\y B {k,M)\ c dv. 



As before, we generate the total power spectra (Prr(A:) and 
PBB(k)) by taking the sum of these parts, 



Prr(*0 = 

m = 



PKK(k) + P£ R (k) 

PBB(k)+P£ B (k). 



(25) 



To generate the cross-power spectrum, we take one fac- 
tor from each of the different power spectrum terms. This 
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Table 1. Fiducial model parameters for power spectra in Figure I 



0.001 0.01 0.1 I 10 100 

k [h Mpc"'] 

Figure 2. Power spectra for red (-PRR(fc)) and blue (PBB(k)) 
galaxies and their cross-power spectrum (-PrbCO) compared to 
the dark matter (Pdd(&)) and galaxy (Pgg(^)) power spectra. 



makes the halo-halo term 



LIN 



(fc) 



_P_ 

P_ 

fin 



fn{v)^-b{v)y^{k,M)dv 



M{v) 



b(v)yB{k,M)dv 



(26) 



For the Poisson term, we simply replace the second moment 
of the galaxy number relations with the product of (N)r and 
(N)b, 



(N)r(N) i 



\y{k,M)\ c dv, 



(27) 



PnB(k) = P fiy) 

J (27r) J n R nM(^) 

where f(u) and y(k,M) are the geometric means of the red 
and blue values. 

The results of such a calculation are shown in Figure ^| 
For our fiducial model, we choose the set of input parame- 
ters in a ACDM model given in Table 1. As we would gen- 
erally expect, the red galaxies show a stronger biasing than 
either the total sample or the blue sample, as well as trac- 
ing the shape of the dark matter power spectrum (-Pdd(&)) 
more closely. The blue galaxies are anti-biased relative to the 
normal galaxy power spectrum, and demonstrate a slightly 
steeper slope. Additionally, the blue galaxies demonstrate a 
sharp break from a power law at small scales. This effect 
is due solely to the (N)(M) relation for the blue galaxies 
and not the halo profiles; the larger number of blue galaxies 
in smaller, less massive halos sets in at this scale, driving 
the power up. Remarkably, however, the galaxy populations 
that generate these power spectra combine to produce a to- 
tal galaxy power spectrum with simple power law behavior. 
The exact comparison of these predicted power spectra to 
those from simulations we leave as a detail for future work; 
for now we are more interested in the flexibility of the model 
than precise values for parameters. 

With this machinery in place, we can calculate the rel- 
ative bias between the various power spectra: 



&rb(A:) 



PrrW 
P BB (fc) 



bnG(k) 



2 _ PrrW 
P GG (k) 



(28) 



Description 


Parameter 


Value 


Red Unit Mass Scale 


Mro 


3 x lO 12 ft- 1 M 


Blue Unit Mass Scale 


M bo 


7 x 10 13 h- 1 MQ 


Red Mass Scaling Index 


7b 


0.9 


Blue Mass Scaling Index 


7b 


0.7 


Gaussian Normalization 


A 


0.5 


Gaussian Mass Scale 


m Bb 


11.75 


Outer Galaxy Ratio 


V 


4 


Inner Galaxy Ratio 


l-l 


10 
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1 1 1 - 
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■ — v« 

\bg« 

J i 1 i 





0.001 0.01 
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Figure 3. Relative biases between red and blue galaxies (6rb(^)) 
and the red, blue and cross power spectra and the galaxy power 
spectrum (6 RG (fe), bBG( fc ): and &RBG( fc ), respectively). 



, ,,s2 fBEi(fc) 

bBG(k) = 



, /,x2 PRB(k) 

bRBG(k) = 



P GG (k) ~" ' P GG (k) 

We choose relative biases between the various galaxy power 
spectra rather than the absolute biases relative to the dark 
matter for two reasons. First, in each of the cases in §[l| 
where we vary a parameter in our model, at least one of 
P GG (k), -pRR(fc) or Pbb(A;) remains roughly constant, so we 
can use that power spectrum as a baseline for seeing how 
the other one or two vary. Second, while the absolute bi- 
ases can be measured using galaxy magnification bias, the 
relative biases involve real clustering that can be measured 
over a much wider range of redshift for a given photometric 
or spectroscopic survey (the evolution of these biases over 
redshift will be left for future work). 

As with the results calculated by Seljak, the relative bi- 
ases shown in Figure ^ are constant on large scales, whereas 
on small scales there is considerable variation, particularly 
in the &RB(fc) and bB G (k) biases. As we will see later, the 
behavior of these biases is a strong function of the model in- 
put, suggesting that reasonably small error bars on the the 
bias in wavenumber space could act as a powerful constraint 
on the model parameters. 
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Figure 4. bjiQ(k) for several values of iirj in the limit where 
A^BO 3> Mro- J'ggW an d PBB(k) are unaffected by changes in 
/it] in this limit. 
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b RG (k), n = 10,ti = 4 

b BG (k),n=10,r,=4 

b RO (k),n = 20,r, = 2 

- b^fk),^ 20,^ = 2 



Ll I I I I I I 

0.001 0.01 0.1 1 10 100 

kfhMpc"'] 

Figure 5. 6r,g(&) and JbgC 1 ) f° r equal values of /irj in the Mro ~ 
Mbo regime. Pgg(^) remains constant for all values of firj by 
construction. 



4 VARIATIONS 

4.1 Inner and Outer Ratios 

As Equations |l^ and suggest, the effect of our choices 
for rj and fj, will be at least partially determined by the ra- 
tio of Mro/Mbo- In the regime where Mbo 3> Mro as m 
the fiducial model, the actual values of fj. and r\ are not so 
important to the resulting power spectra as their product. 
Likewise, since the profile that the blue galaxies will popu- 
late is relatively flat in the center of each halo, the effects 
of changing the value of /irj does not significantly change 
the clustering of the blue galaxies; the majority of the mass 
associated with normal NFW profiles is outside the p ~ r a 
region anyway and, since lvb > ct by construction, this will 
be even more pronounced for the blue galaxies. Thus, we do 
not expect PBB(fc) to vary significantly with \ir\ and we have 
constructed or and cvb such that Pgg(&) will not vary, so 
the only sensitivity to prj we expect to see is in PRR(fc) and 
ft.B(fe) and the associated relative biases. In Figure 0, we 
show 6rg(/c) for several different values of fir). There is some 
change in the shape of the bias, perhaps indicative of a more 
negative qr in the high fir) models leading to a greater pop- 
ulation of the inner regions of the halos with red galaxies. 
Indeed, in most of these models, the matter associated with 
red galaxies only exceeds that associated with blue galaxies 
in the very inner regions of the halo. 

In a model where Mbo ~ Mro, the degeneracy between 
fi and i) is broken. Figure U shows 6rg(&) and 6bg(&) for two 
models where Mbo = M R0 = 3 x 10 12 h' 1 M Q and fit) = 40. 
Since Equations |l^ and [22] are no longer dominated by the 
ratio of AIro/Mbo, we can see significant shifts in the biases 
of both red and blue galaxies relative to Pgg(&)- In this 
case, we have a more equal distribution between the mass 
associated with red and blue galaxies and less extreme rel- 
ative halo profile normalizations. Clearly, using future mea- 
surements to constrain these parameters will require using 
multiple biases to minimize these degeneracies. 



4.2 (N)(M) Relations 

Unlike modifying the halo profiles, changing the parameters 
in the (N)(M) relations can have significant effects on the 
shapes of all the power spectra not just Prr(/c) and PBs(k). 
The general effect of each of the modifications is to change 
the behavior of the Poisson term in the power spectra, but 
the specific effects for each modification show considerable 
and sometimes surprising variations. 

We begin with the unit mass scale for the blue galax- 
ies, Mbo- In general, this parameter does not strongly af- 
fect the total galaxy power spectrum; there is some slight 
variation (~ 5%) over the range 10 12 /i _1 Mq < Also < 
W 14 h~ 1 Mq in the quasi- linear regime of the power spec- 
trum (k ~ hMpc -1 ). However, there is significant change 
in the relative biases of the red and blue galaxies, as shown 
in Figure ^. Additionally, as log (Mbo) approaches Mb s , the 
non-power law behavior of PBB(fc) is considerably damped 
out. 

In contrast to A^bo, modifying the values of Mro leads 
to large variations in the shape and amplitude of Pgg(&), 
with lesser amplitude shifts to Prr(&) and almost no effect 
on PBB(fc) (as we might suspect). Figure (?] shows &RB(fc) 
and &bg(&) for a two decade range in Mro values. As Mro 
approaches A/bo, the effect of the Gaussian term in (N)b on 
Pgg(^) increases, leading to a galaxy power spectrum with 
a strong break in its power law at large wavenumber. Con- 
versely, at lower Mro, the galaxy power spectrum inflects, 
leading to a stronger anti-bias in PBB(fc) relative to Pgg(&) 
around k ~ 1. Additionally, in this limit we can see the ef- 
fect of the Gaussian component in the transformation from 
77 to 77' (Equation [il] ) changing the effective value of fir/. 

Due to the sub-dominant role of (N)b to that of (N)r 
over most of the mass range, varying 7b does not signifi- 
cantly change any of the power spectra. Changing 7r alters 
PRR.(fc) and Pgg(&) slightly, generally smoothing out the 
variations in &RG(fc) over k for larger values of 7r. 

Changing A and Mb s has minimal effect on the red 
power spectrum, so long as we are in the Mro, Mbo 2> 10 Mbb 
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regime. There are, however, effects on Pgg(&) and PBB(k) 
(quite strong effects in the case of Mb s ) and we can begin 
to see our approximation neglecting the Gaussian term in 
Equation ^ break down as Mb s approaches M» , resulting in 
biased Pbb (k) relative to Pgg (k) at small scales. Increasing 
the value of A has a slight global effect on &rb (k) and 6rg (k) 
(Figure but mostly it controls the on-set of the break 
in the Pbb(&) power-law with large values of A lead to a 
sharper break. The mass scale for the contribution of the 
Gaussian term in (N)b plays a much more significant role. 
Large values of Mb s increase PBB(fc) (and, to a lesser extent, 
Pgg(&)) on all scales, leading to a suppression of 6rb(&) and 
&RG(fc) as Mbb increases (Figure ^). Likewise, as the mass 
of the halos experiencing this boost in P£&(k) increases, the 
onset of the bump in the power law for PBB(fc) occurs on 
larger and larger scales. 



As mentioned above, the over-riding theme of these vari- 
ations appears to be the effect of changing at what scale and 
in what way the Poisson contribution to the power spectra 
sets in. As Mbo and Mro increase, fewer galaxies resulting 
from the power law parts of (N) B and (N) R populate the 
lower mass halos and small scale power-law behavior is sup- 
pressed in favor of the Gaussian contribution to (N)b. The 
degree to which this non-power law behavior in the red and 
blue galaxies is recreated in the data and simulations should 
be an excellent clue as to setting the relative amplitude of 
Mbo and Mro as well as Mb s - 



5 CONCLUSION 

We have shown that, through relatively simple modifications 
to the (N) and halo profile relations in the standard formal- 
ism, we can generate reasonable power spectra for red and 
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blue galaxies, as well as a number of relative biases between 
the various galaxy power spectra. By manipulation of the 
parameters constituting the (N) relations, we have consid- 
erable ability to modify the shapes of the red and blue power 
spectra. Likewise, in the limit where we are not dominated 
by the (N) relations, our choice of halo profile parameters 
allows us to set the relative large-scale biases between the 
various power spectra over a large range, while keeping the 
shapes of the power spectra relatively constant. There is 
some degree of degeneracy between the choice of Mbo & Mro 
and the halo profiles for each of the galaxy sub-populations, 
suggesting that a simple fit to (N) from simulations without 
including this effect could result in an apparently larger dif- 
ference between the two mass scales than the data actually 
indicate. Still, the degree of flexibility within the framework 
and sensitivity to the various input parameters suggest that 
measurements of these relative biases using current large 
galaxy survey data will provide strong constraints on these 
input parameters and the outputs of simulations. 
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